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I .  ABSTRACT 


The  purpose  of  this  project  was  to  perform  theoretical  and 
experimental  investigation  on  reactive  and  nonreactive  collision 
processes  that  occur  in  chemical  laser  systems.  The  theoretical 
investigation  was  focused  on  developing  reliable  and  tractable 
numerical  techniques  that  describe  the  dynamics  of  the  collision 
processes.  The  experimental  effort  was  devoted  to  a  chemical  laser 
scheme  based  on  gas  phase  chemiluminescence  reactions,  which  allows 
premixing  of  the  reagents. 

II.  SCIENTIFIC  BACKGROUND 

Chemical  lasers  (for  example,  HF,  DF,  and  iodine)  have 
efficiently  achieved  power  outputs  of  hundreds  of  kilowatts.  So 
far,  our  understanding  of  these  lasers  is  hampered  by  our  inability 
to  predict  and  interpret  the  energy  distribution  in  the  collision 
processes  that  pump  the  lasing  states.  The  difficulty  is  not  so 
much  in  understanding  the  physical  laws  that  govern  these  processes 
as  in  solving  the  equations  that  exactly  describe  them;  the  equations 
are  much  too  complicated  to  be  solvable.  Therefore,  researchers 
working  with  chemical  laser  systems  would  like  to  see  the  development 
of  a  technique  to  solve  the  equations  that  is  both  accurate  and 
numerically  tractable.  Since  the  exact  quantum-mechanical  approach 
is  not  numerically  tractable  in  general,  one  must  rely  on  various 
approximate  results.  The  most  popular  approximate  technique  for 
describing  molecular  collisions  is  the  classical  trajectory  method, 
in  which  the  dynamics  of  collisions  are  described  completely  by 


classical  mechanics  rather  than  quantum  mechanics.  When  this 
technique  is  combined  with  the  commonly  used  procedure  for  classify¬ 
ing  trajectories  as  to  their  initial  and  final  states,  it  is  often 
referred  to  as  the  quasi-classi cal  (QC)  method.  Although  the  QC 
method  is  easy  to  apply  and  interpret,  it  is  often  not  sufficiently 
accurate  in  describing  microscopic  state-to-state  aspects  of  colli¬ 
sions  which  are  important  in  the  study  of  chemical  laser  systems. 
Furthermore,  because  the  method  is  intrinsically  classical,  it  is 
not  capable  of  describing  classically  forbidden  processes  resulting 
from  quantum  tunneling.  To  overcome  this  weakness  of  the  QC  method, 
a  method  usually  called  a  semiclassical  (SC)  method  was  introduced. 

In  this  method,  classical  trajectories  are  allowed  to  take  on  com¬ 
plex  values,  thus  enabling  one  to  approximately  describe  classically 
forbidden  processes.  However,  one  here  carries  the  burden  of  working 
with  complex-valued  trajectories,  which  are  very  difficult  to  handle. 
It  therefore  seems  highly  desirable  to  develop  a  new  technique  that 
combines  to  a  maximum  degree  the  simplicity  of  the  QC  method  with 
the  accuracy  of  the  fully  quantum-mechanical  (QM)  method.  This  is 
where  our  theoretical  effort  under  Grant  No.  AFOSR-79-0122  has  been 
directed. 

Chemical  lasers  such  as  111',  HC1 ,  anil  CO  operate  on  vibrational 
or  rotational  transitions  and  arc  thus  infrared.  A  recent  active 
field  of  research  has  been  the  development  of  visible  wavelength 
chemical  lasers  that  operate  on  electronic  transitions.  To  achieve 
lasing  on  an  electronic  transition  in  a  typical  diatomic  molecule, 
it  is  necessary  to  rapidly  generate  a  supply  of  molecules  in  the  upper 
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level  of  the  laser  transition.  At  visible  wavelengths  this  requirement 
is  difficult  to  satisfy  with  a  traditional  mixing-type  chemical  laser 
that  employs  chemical  reactions  which  initiate  spontaneously  upon 
mixing  of  the  reagents.  An  alternative  scheme  that  we  have  pursued 
involves  gas-phase  chemiluminescence  reactions  that  do  not  initiate 
until  "triggered"  ’  .  a  suitable  low-power  ignition  mechanism.  Thus, 
reagents  can  be  premixed,  and  a  small  input  of  optical  or  electrical 
energy  can  be  used  to  produce  an  initial  chain  carrier  species.  The 
resulting  chemical  reaction  is  thus  self-sustaining  via  branched 
chain  reaction  kinetics.  Our  experimental  effort  was  devoted  to 
this  scheme,  utilizing  chemiluminescence  reactions  that  are  triggered 
by  a  pulsed  laser. 

III.  OVERVIEW  OF  THE  WORK  PERFORMED 

In  research  supported  by  Grant  AFOSR-79-0122 ,  we  have  accom¬ 
plished  the  following: 

1.  Developed  a  new  approach  to  molecular  collisions  that 
appears  to  have  overcome  several  difficulties  of  the  QC  and  SC 
methods . 

2.  Modified  the  SC  method  so  as  to  accurately  describe  near- 
rcsonancc  vibrat ional-to-elcctronic  energy  transfer  processes. 

3.  Obtained  intense  blue-green  chemiluminescence  from  the 
d*TT  -►  a3ir  transition  of  C?  by  premixing  CIO?  with  acetylene  (C?ll?)  and 
triggering  t he  reaction  with  an  XoF  excimer  laser. 

Progress  in  the  above  three  areas  is  described  in  the  follow¬ 
ing  subsections.  Further  details  can  be  found  in  the  appended  reports. 
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A.  A  New  Approach  to  Molecular  Collisions:  Statistical 
Quasi-Classical  Method 

Our  theoretical  effort  was  focused  on  developing  a  technique 
that  combines  to  a  maximum  degree  the  simplicity  of  the  QC  method 
and  the  accuracy  of  the  QM  methou.  This  of  course  is  the  spirit 
with  which  the  SC  method  was  developed.  To  avoid  complications  due 
to  the  handling  of  complex-valued  trajectories,  we  elected  to  keep 
the  trajectories  real,  thus  describing  the  dynamics  of  collisions 
completely  by  classical  mechanics.  However,  we  used  a  quantum 
statistical  representation  of  initial  and  final  states  instead  of  the 
usual  representation  based  on  the  classical  probability  distribution 
that  is  used  in  the  QC  method.  Thus,  our  new  method  may  be  con¬ 
sidered  to  be  the  classical  trajectory  method  with  a  touch  of  quantum 
statistics  in  representing  initial  and  final  states.  For  this  reason, 
we  call  our  method  the  statistical  quasi-classical  (SQC)  method. 

The  quantum  statistical  representation  of  the  states  seems  to 
be  most  conveniently  achieved  by  using  the  Wigner  distribution  func¬ 
tion, 

W(q,p,t)  =  —-^dx  expC-2ipx/h)il)*(q-x,t)\Kq  +  x,t)  .  (1) 

Roughly  speaking,  the  Wigner  function  gives  the  quantum-mechanical 
probability  distribution  at  time  t  in  phase  space  (q,p)  for  the 
system  described  by  the  wavcfunction  i{i(q,t}.  To  use  it  in  our  SQC 
method,  we  first  need  to  find  the  initial  Wigner  function 
W(q,p,t  =  -">)  corresponding  to  the  initial  state  <|i(q,t  -  -cc’) .  Fateh 
element  (q,p)  of  the  initial  Wigner  function  is  then  propagated  in 
time  using  the  classical  Hamilton  equations  of  motion.  This 


5 


propagation  is  achieved  in  practice  by  representing  the  initial  Wigncr 
distribution  by  a  finite  number,  N,  of  points  on  an  appropriate  two- 
dimensional  grid  in  the  (q,p)  plane  and  then  integrating  the  Hamilton 
equations  of  motion  for  each  of  the  N  points  (q,p)  on  the  grid.  The 
collection  of  the  final  values  of  q  and  p  obtained  after  the  integra¬ 
tion  constitutes  the  final  Wigner  distribution  function  W(q,p,t  =  «) , 
from  which  the  transition  probability  can  be  obtained. 

Since  the  SQC  method  contains  some  quantum-mechanical  features 
not  contained  in  the  QC  method,  one  might  expect  it  to  be  able  to  at 
least  partially  describe  classically  forbidden  processes.  The  chief 
disadvantage  of  the  SQC  method,  on  the  other  hand,  is  that  it  requires 
a  large  number  of  trajectories  to  be  computed,  as  each  element  (q,p) 
of  the  Wigner  distribution  must  be  propagated.  However,  this  is  only 
a  small  disadvantage,  considering  the  advantage  of  not  having  to  deal 
with  complex- valued  trajectories. 

The  SQC  method  was  applied  to  the  collinear  inelastic 
collision. 

He  +  H2  Cv^ )  -*■  He  +  H2  (v^O  .  (2) 

Calculations  clearly  showed  that  the  SQC  method  is  more  accurate 
than  the  QC  method.  They  also  indicated  that  the  SQC  method  is, 
or  at  least  has  a  potential  to  be,  capable  of  accurately  describing 
classically  forbidden  processes.  It  appears,  therefore,  that  the  SQC 
method  represents  an  appealing  alternative  to  the  existing  tech¬ 
niques  for  describing  molecular  collisions,  especially  when  more 
accurate  information  is  desired  than  can  be  obtained  by  the  QC  method 
and  when  the  fully  QM  method  or  the  SC  method  is  too  difficult  to 
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handle. 


B.  Modification  of  the  Semiclassical  Method:  Trajectory-Based 
Decoupling  Scheme 

The  major  difficulty  of  the  SC  method  in  describing  collision 
processes  lies  in  the  complexity  of  integrating  complex-valued  tra¬ 
jectories.  One  of  the  approximate  schemes  designed  to  alleviate 
this  difficulty  is  the  "decoupling"  scheme,  in  which  trajectories 
are  assumed  to  take  on  real  values  except  when  they  hit  (avoided) 
crossing  points  at  which  electronic  transitions  may  occur  in  the  SC 
theory.  The  scheme  then  assumes  that  all  but  one  appropriately 

chosen  coordinate  n  and  momentum  P  can  be  "frozen"  at  the  crossing 

n  b 

points  and  calculates  the  transition  probability  in  a  Stuckelberg 
manner,  allowing  only  the  "decoupled"  coordinate  n  and  momentum 
to  take  on  complex  values.  It  is  important  to  make  an  intelligent 
choice  of  the  decoupled  coordinate  n  if  the  approximation  is  to  be 
reliable.  The  usual  way  of  choosing  it  is  to  look  at  the  structure 
of  the  potential  energy  surface  in  the  neighborhood  of  choosing  points 
and  determine  which  component  of  momentum  is  most  effective  in  leading 
to  electronic  transitions.  For  example,  for  the  colinear  collision 

X*(2P^  +  H2(v.)  -  X(2p3/2)  +  H2(vf),  (3) 

where  X  is  cither  F  or  Br,  the  decoupled  coordinate  n  was  chosen  to 
be  the  translational  coordinate  R.  Although  SC  calculations  with  this 
choice  of  n  yielded  energy  transfer  probabilities  in  good  agreement 
with  the  corresponding  QM  calculations  for  X  =  F,  a  strong  disagree¬ 
ment  was  found  for  X  =  Br.  This  was  attributed  to  the  inability  of 
the  SC  method  to  describe  near-resonance  processes  involving 
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vibrational  or  rotational  energy  levels,  because  internal  (vibra¬ 
tional  and  rotational)  degrees  of  freedom  are  treated  classically 
in  the  SC  method.  However,  it  is  not  entirely  clear  that  resonance 
effects  should  be  absent  simply  because  internal  motions  are  treated 
classically.  A  good  agreement  might  be  obtained  with  a  better  choice 
of  the  decoupled  coordinate,  even  if  we  use  the  SC  method  or  even 
the  QC  method.  This  is  an  important  consideration  because  it 
determines  whether  one  can  use  a  numerically  tractable  technique  (the 
QC  method)  to  describe  near-resonance  processes  of  importance  in 
chemical  laser  systems. 

A  better  choice  of  the  decoupled  coordinate  n,  we  believe, 
can  be  achieved  by  taking  n  to  be  the  coordinate  along  the  direction 
of  motion  at  a  crossing  point.  The  reason  for  this  particular  choice 
can  be  found  in  the  reprint  of  Appendix  B  (Lee,  I.am,  DcV'ries,  and 
George) .  Here  we  only  note  that  n  determined  by  the  above  choice 
is  generally  different  for  different  trajectories,  and  therefore  has 
to  be  defined  for  each  trajectory  at  each  encounter  with  a  crossing 
point.  For  this  reason,  the  new  decoupling  scheme  is  called  the 
trajectory-based  decoupling  scheme.  QC  calculations  have  been  car¬ 
ried  out  for  the  near-resonance  vibrational-to-elcctronic  energy 
transfer  process 

Br(2P3/2)  +  H2 Cvj  =  1)  -  Br(2P,^)  +  H?(vf  =  0)  (4) 

using  the  trajectory-based  decoupling  scheme.  The  agreement  with  the 
QM  calculations  is  much  better  witli  the  trajectory-based  scheme  than 
with  the  old  scheme.  Furthermore,  as  we  artificially  varied  the  spin- 
orbit  splitting  of  Br  around  the  fixed  vibrational  frequency  of  II^, 
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the  new  scheme  clearly  reproduced  the  resonance  peak  structure  of 
the  transition  probability,  whereas  the  old  scheme  failed.  It 
therefore  seems  that  our  new  scheme  opens  the  door  for  the  QC 
method  to  describe  near-resonance  energy  transfer  processes. 

The  work  described  in  this  section  was  carried  out  in  co¬ 
operation  with  Prof.  Thomas  F.  George  and  his  coworkers  at  the 
University  of  Rochester. 

C.  Laser- Initiated  Chemiluminescence  Experiments 

To  obtain  chemiluminescence  from  C2  using  the  premixing 
technique,  we  first  generated  chlorine  dioxide  (C102)  by  passing 
Cl2  gas  over  NaC102  crystals.  The  C102  was  then  collected  into  the 
reaction  cell,  to  which  acetylene  (C21I2)  was  added  to  allow  it  to 

be  mixed  with  the  C102.  A  single  10-mJ,  20-nsec  pulse  of  3500  A 

light  from  a  XeF  cxcimer  laser  was  then  used  to  detonate  the  mixture 
The  detonation  resulted  in  an  intense  blue-green  chemiluminescence, 
which  was  easily  identified  as  the  C  2  ( d  3  m  ->•  a3 tt )  Swan  bands.  From 
the  observation  of  the  temporal  behavior  of  the  chemiluminescence, 
we  found  that  the  rise  time  of  the  pulse  is  typically  30  to 

50  nsec,  while  the  "burn-time"  exceeds  1  msec.  The  time  delay  befor 

the  rapid  onset  of  light  emission  was  on  the  order  of  the  time  re¬ 
quired  for  a  shock  front  to  travel  from  the  laser  input  window  to 
the  side-arm  where  light  emission  was  observed.  We  do  not  know 
the  exact  mechanism  of  reaction  and  therefore  cannot  speculate  in 
detail  regarding  the  kinetics  of  the  reaction.  The  blue-green 
chemiluminescence  from  t lie  C2  Swan  band  was  .also  observed  from 
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HN3  +  C2II2  +N20  +  ile  mixtures  using  a  single  20-mJ,  20-nsec  pulse  of 
2490-X  KrF  excimer  laser  light.  In  addition,  strong  CN  violet  band 
emission  was  also  observed  from  this  mixture.  To  determine  whether 
a  chemically  pumped  C2  laser  is  possible  from  the  schemes  described 
above,  however,  further  investigation  is  necessary  on  the  kinetics  of 
the  reactions  involved. 

The  work  described  in  this  section  was  carried  out  in  coopera¬ 
tion  with  Dr.  David  Benard  at  AFWL. 

IV.  PUBLICATIONS 

The  following  publications,  produced  in  conjunction  with  this 
work,  are  reproduced  as  appendixes  to  this  report: 


Hai-Koong  Lee  and  Marian  0.  Scully,  "A  new  approach  to  molecular 
collisions:  Statistical  quasiclassical  method,"  J.  Chem.  Phvs . 
73(5):  223S-2242 ,  1  Sept.  1980. 

Hai-Koong  Lee,  Kai  S.  Lam,  Paul  L.  DeVries,  and  Thomas  F.  George, 

"A  new  semiclassical  decoupling  scheme  for  electronic  transitions 
in  molecular  collisions:  Application  to  vibrational-to-clectronic 
energy  transfer,"  J.  Chem.  Phys .  73(1):  206-212,  1  July  1980. 

M.  D.  Burrows  and  D.  J.  Benard,  "Pulsed  visible  chemiluminescence 
experiments,"  Laser  Digest,  summer  1980. 
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APPENDIX  A 


A  new  approach  to  molecular  collisions:  Statistical 
quasiclassical  methoda) 

Hai-Woong  Lee  and  Marian  O.  Scully 

Department  of  Physics  and  Optical  Sciences  Center,  University  of  Arizona.  Tucson.  Arizona  8572 1 

and  Projek tgru ppe  fur  Laserforschung  der  Max-Ptanck  Gesellschaft,  D-8046  Carefling  bei  Munchen. 

West  Germany 

(Received  25  March  1980;  accepted  12  May  1980) 

A  new  approach  is  presented  in  which  classical  mechanics  is  combined  with  quantum  statistics  to  describe 
molecular  collisions.  In  this  approach,  the  dynamics  of  collisions  is  described  by  classical  trajectories  as  in  the 
widely  used  quasiclassical  method.  However,  initial  and  final  internal  states  are  represented  in  phase  space  in 
a  quantum  statistical  way,  using  the  Wigner  distribution  function.  Results  of  calculations  performed  on  a 
collinear  He-H.  collision  indicate  that  this  new  method  is  more  accurate  than  the  quasiclassical  method, 
especially  when  the  initial  vibrational  energy  is  low.  Moreover,  the  new  method  is  capable  of  describing 
classically  forbidden  processes  that  cannot  be  accounted  for  by  the  quasiclassical  method. 


I.  INTRODUCTION 

The  classical  trajectory  method  is  probably  the  most 
popular  technique  for  describing  the  dynamics  of  molec¬ 
ular  collisions. 1,2  As  the  name  of  the  method  implies, 
the  colliding  molecules  obey  Hamilton’s  equations  of 
motion.  The  technique  is  often  referred  to  as  the  quasi¬ 
classical  (QC)  method  when  combined  with  the  common¬ 
ly  used  procedure  for  classifying  trajectories  as  to  their 
initial  and  final  states.  The  strength  of  the  QC  method 
is  that  it  is  easy  to  apply  and  requires  less  computer 
time  and  expense  than  the  accurate  quantum -mechanical 
(QM)  method.  However,  the  QC  method  is  often  not  suf¬ 
ficiently  accurate  in  describing  microscopic  state-to- 
state  aspects  of  collisions.  Furthermore,  the  QC  meth¬ 
od,  due  to  its  intrinsically  classical  nature,  is  not  ca¬ 
pable  of  describing  classically  forbidden  processes  re¬ 
sulting  from  quantum  mechanical  tunneling.  Therefore, 
it  is  highly  desirable  to  develop  a  technique  that  over¬ 
comes  the  weakness  of  the  QC  method,  while  retaining 
its  strong  points  as  much  as  possible.  Progress  has 
been  made  along  this  direction  in  recent  years.3-5 

In  this  paper,  we  introduce  an  approach  termed  the 
statistical  quasiclassical  (SQC)  method,  which  is  essen¬ 
tially  the  classical  trajectory  method  combined  with  a 
quantum  statistical  representation  of  initial  and  final 
states.  The  idea  of  the  method  stems  from  the  fact  that 
the  usual  assignment  of  initial  conditions  in  the  QC  meth¬ 
od  is  often  inaccurate  from  the  quantum  mechanical 
point  of  view.  For  example,  the  selection  of  the  vibra¬ 
tional  phase  of  a  harmonic  oscillator  is  based  on  the 
probability  distribution  resulting  from  the  classical  har¬ 
monic  motion.  This,  however,  approaches  the  correct 
quantum  mechanical  distribution  only  in  the  limit  of  high 
vibrational  energy.  One  might  expect  to  obtain  higher 
accuracy  by  using  an  initial  assignment  based  on  the 
quantum  mechanical  nature  of  the  system  being  investi¬ 
gated.  This  consideration  immediately  leads  to  the  idea 
of  using  an  appropriate  probability  distribution  function 
In  phase-space,  derived  using  quantum  mechanical  wave 
functions.  We  choose  for  this  purpose  the  Wigner  dis- 
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tribution  function,  ®*7  which  has  been  widely  used  in  the 
past  for  describing  a  variety  of  systems. a-u  Thus,  in 
our  SQC  method,  initial  internal  states  are  represented 
by  the  corresponding  Wigner  distribution,  which  then  is 
propagated  in  time  using  the  classical  trajectory  method. 
The  resulting  final  Wigner  distribution  carries  informa¬ 
tion  on  the  final  state,  and  therefore  can  be  used  to  ob¬ 
tain  transition  probabilities.  The  Wigner  distribution 
function  yields  the  correct  quantum  mechanical  proba¬ 
bility  distribution  for  both  position  and  momentum  of  the 
system  under  consideration.  Therefore,  the  SQC  method 
contains  at  least  some  quantum  mechanical  features  not 
contained  in  the  QC  method.  In  this  regard,  it  is  of  pri¬ 
mary  interest  to  see  how  well  the  SQC  method  can  de¬ 
scribe  classically  forbidden  processes.  In  Sec.  II  we 
give  a  brief  review  of  the  basic  properties  of  the  Wigner 
distribution  function,  useful  for  the  development  of  the 
SQC  method.  Section  III  follows  with  a  description  of  the 
SQC  method,  and  its  application  to  a  collinear  He-H2 
collision  is  described  in  Sec.  IV. 

II.  THE  WIGNER  DISTRIBUTION  FUNCTION 

The  Wigner  distribution  function  W(.q,p,t )  may  be  de¬ 
fined  as  the  Fourier  inverse  of  the  characteristic  func¬ 
tion 

C(Z,T],t)=  f  dq  (/>*  (q ,  t)  exp[i(£?  +  T)p)]i>(q ,  t) 

=  Tr{p(q,/>,f)exp[iU5  +  7jp)]}  ,  (1) 

i.e., 

W(q,p,t)  =  f  dt,  j  drj  exp[-  i(£<7  +  »)/>)] 

X  ^ dq  (q ,t)  exp[i(l,q  +qp)]il’(q ,  t)^ 

=  Tr[p(q,p,i)  6{q -q)  6(p -/>']  .  (2) 

Here  £  and  r\  are  real  parameters,  q  and  p  denote  coor¬ 
dinate  and  momentum  operators,  and  q  and  p  corre¬ 
sponding  eigenvalues.  It  can  be  easily  shown7  that  Eq. 

>2)  reduces  to 

W(<J  ,/>,()  =  -“  J  d\  {q  -  x  ,t)  iMq  u.I)  ,  (3) 
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an  expression  first  given  by  Wigner.9  If  we  expand  $  in 
terms  of  an  orthonormal  set  of  eigenfunctions 

=  X  <>„(')*.(?)  ,  (4) 

IV 

Eq.  (3)  becomes 


a  * 

where 

p)  dxex p(~^jf£)  -  *>  +  *)•  (6) 

It  can  also  be  shown1  that  the  functions  uin „’s  form  a  com¬ 
plete  orthogonal  set  in  the  phase  space  (q,p)  , 

jdq  j dp  w'jq ,  p)  u-V, .(</,/>)  =  >  <7a) 


-q'q~p~p,)  ■  (7b) 


If  we  integrate  the  function  il  given  by  Eq.  (3)  over 
p,  we  obtain 


fdpW(q,p,t)=4*(q,t)t(q  ,t) 


(8a) 


Similarly,  we  get 

JdpW(q,p,l)=i*(p,t)iip,t)  ,  (8b) 

where  ?(/),/)  is  the  Fourier  inverse  of  it(q,t),  or  the 
wave  function  in  the  momentum  space.  Equation  (8)  in¬ 
dicates  that  the  Wigner  distribution  function  yields  the 
correct  quantum  mechanical  probability  distribution  in 
both  coordinate  and  momentum  spaces. 


III.  STATISTICAL  QUASICLASSICAL  METHOD 

In  this  section,  we  describe  our  new  approach,  the 
SQC  method.  In  principle,  this  method  can  be  applied 
to  any  type  of  collision  process.  For  convenience  of 
presentation,  however,  we  consider  a  nonreactive  atom- 
diatom  collision 

A  +  BC(n)—  A  +  BC(m)  ,  (9) 

where  n  and  m  refer  to  any  relevant  set  of  initial  and 
final  quantum  numbers.  We  denote  the  translational 
coordinates  and  momenta  collectively  by  Q  and  P  and 
the  internal  coordinates  and  momenta  by  q  and  p.  At 
the  initial  time  t  =  -•*,  the  diatom  DC  is  prepared  in  its 
nth  eigenstate 

tack,  f  =  -°°)=  </’„(?>  •  (10) 

In  our  approach,  this  initial  state  is  represented  by  the 
corresponding  Wigner  distribution  function,  which  can 
be  written  by  using  Eqs.  (5)  and  (6)  as 

»'(<?, P,  /  =  -«>)  =<»•,,„(</,/>)  -  ~  j  dx cxp(— 

XC(<7  -*)<i’„(<7  •  v)  .  (11) 

The  initial  translational  coordinates  and  momenta  are 
specified  as  usual:  the  initial  Q  corresponds  to  an  arbi¬ 
trarily  large  separation  of  A  apd  DC  at  which  the  A -DC 


interaction  is  negligible,  and  the  initial  P  is  determined 
by  the  collision  energy.  Now,  our  system,  A  +  DC,  ini¬ 
tially  specified  by  Q,  P  and  W(q,p,  t  =  -°°),  evolves  in 
time  according  to  Hamilton’s  equations  of  motion.  To 
integrate  Hamilton’s  equations,  we  first  represent  the 
initial  distribution  W{q  ,p,  /  =  -*>)  by  a  finite  number  (N) 
of  points  on  an  appropriate  two-dimensional  grid  in  the 
iq,p)  plane.  The  integration  is  then  carried  out  N  times 
for  fixed  initial  values  Q,  P  and  for  each  point  (q,p)  on 
the  grid.  The  collection  of  final  values  of  q  and  p  ob¬ 
tained  after  the  integration  constitutes  the  final  distribu¬ 
tion  function 

W(q,p,  f  =  ®)  =  X)  X  «?<*=•)  *iM'p)  .  (12) 

<  l 

corresponding  to  the  final  state  of  the  diatom  BC, 

/  =«>)  =XT  "  =°°)<Mv)  •  (13) 

i 

The  transition  probability  «—  w  can  then  be  obtained  from 

P»-,  =  l'i-«=,0)|*  =  2ir«  jdq  j dpumJq,p)W(q,p,  t  =»)  , 

(14) 

where  the  orthogonality  relation,  Eq.  (7a),  has  been 
used  in  obtaining  Eq.  (14).  Thus,  the  recipe  to  obtain 
the  transition  probability  P„.n  using  the  present  SQC 
method  is  simple:  one  first  finds  the  initial  Wigner  dis¬ 
tribution  from  Eq.  (11)  and  then  obtains  the  final 
Wigner  distribution  W{q,p,  I  =*>)  by  integrating. Hamil¬ 
ton’s  equations  of  motion  for  each  element  (q,p)  of 
The  transition  probability  can  then  be  obtained  by  Eq. 

(14). 

It  should  be  noted  that  the  statistical  representation 
of  the  initial  state  according  to  Eq.  (11)  is  what  differ¬ 
entiates  the  SQC  method  from  the  usual  QC  method. 

The  uncertainty  relation  required  by  quantum  mechani¬ 
cal  considerations  is  contained  in  this  statistical  repre¬ 
sentation,  and  for  this  reason  the  SQC  method  is  ex¬ 
pected  to  be  able  to  describe  classically  forbidden  pro¬ 
cesses.  It  should  also  be  noted  that  in  our  method  the 
Wigner  distribution  function  is  used  only  to  represent 
the  initial  and  final  states,  and  the  collision  dynamics 
occurring  between  the  initial  and  final  times  is  still 
governed  completely  by  classical  mechanics.  This  dif¬ 
fers  from  other  works  utilizing  the  Wigner  distribution 
function, 9-11  in  which  the  dynamics  is  governed  by  the 
Schrodinger  equation  of  motion.  One  difficulty  in  usine 
the  essentially  quantum  mechanical  Wigner  distribution 
function  to  determine  initial  and  final  conditions  lor 
purely  classical  equations  of  motion  is  that  some  quan¬ 
tum  mechanical  features  of  the  Wigner  distribution  func¬ 
tion  are  tost  during  the  classical  dynamical  evolution. 
This  is  a  fundamental  limitation  on  the  SQC  method  in 
its  present  form.  This  point  will  be  further  illustrated 
in  the  next  section.  In  spite  of  this  difficulty,  our  de¬ 
sire  to  keep  the  calculational  procedure  at  an  elemen¬ 
tary  level  provides  strong  incentive  for  working  with  the 
classical  equations  of  motion. 

IV.  EXAMPLE 

To  illustrate  our  approach,  we  model  the  colliucar 
He-Hj  collision  as  follows:  ll4  is  approximated  by  a 
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harmonic  oscillator,  and  the  He-H2  interaction  is  an 
exponential  repulsion  between  He  and  the  closer  end  of 
H2.  This  system  has  been  extensively  studied  in  the 
past,  quantum  mechanically13’14  as  well  as  semielassi- 
cally .  li-1‘  We  choose  the  potential  parameters  accord¬ 
ing  to  Secrest  and  Johnson, 13  and  one  should  consult 
their  paper  for  details  of  the  system.  Here  we  only  need 
to  discuss  those  features  that  are  newly  encountered  in 
our  approach. 

Let  us  assume  that  the  harmonic  oscillator  (Hs)  is  ini¬ 
tially  in  its  nth  eigenstate 

ipJq)=NmHn(aq)exp(-2  azq2)  ,  (15) 


where  a  =  ((iiii1/t)1/!;  p  and  u>0  are,  respectively,  the 
mass  and  the  angular  frequency  of  the  harmonic  oscilla¬ 
tor;  JV„  is  the  normalization  constant;  and  Hn  denotes  the 
nth  Hermite  polynomial.  Then  the  initial  Wigner  distri¬ 
bution  function  is  given  by  Eqs.  (5)  and  (6)  as 


»’(?./>.  -°°)  =  un,(q,p)  =  —■  f  dxexp 

xexp[-jo3(<7  -.v)]*/f„(a<7  +  a.v)exp(-5  as(<7  +  x)J]  . 


By  straightforward  algebra,  we  obtain 

exp  [-  (a  v  *  s?f)]  ( 2a  v + 


(16) 


(17) 

or 

wmW.p)  =  ww(,v)=[(-l)"'”('lexp(-y)L.(2y)  ,  (18) 

where  the  energy  parameter  y  is  defined  as 

y  =  azq*  +  pz'aznz  =  {2/hu0){ipu:lqi  +  t>z/2p)  ,  (19) 

and  L„  refers  to  the  nth  Laguerre  polynomial.  For  ex¬ 
ample, 

<*>oo(y)=(l,irfi)exp(-y)  ,  (20a) 

wn(y)  =  (-  1/vH)  (1  -  2y)  exp(-  v)  .  (20b) 


We  note  that  the  initial  distribution  depends  only  on  the 
energy  parameter  y.  This,  however,  is  due  to  the  as¬ 
sumption  that  the  system  is  initially  in  an  energy  eigen¬ 
state.  If  the  system  is  not  in  an  energy  eigenstate,  the 
Wigner  distribution  function  would  contain  off-diagonal 
elements  wu,  i *  j,  which  in  general  depend  on  q  and  p 
separately. 


In  order  to  numerically  represent  the  initial  distribu¬ 
tion  given  by  Eq.  (17),  we  construct  a  two-dimensional 
rectangular  grid  in  the  (q,  p)  plane.  The  size  and  the 
density  of  the  grid  should  be  taken  as  the  smallest  ones 
with  which  the  desired  accuracy  can  lie  achieved.  For 
our  calculations  of  the  He— IL  collision,  for  example, 

20x  20  points  in  the  rectangle  connecting  the  lour  points 
(i  2  'a,  j  2o/i)  gives  sufficiently  accurate  results  for 
n  -  1.  Each  point  (r/(,  />,)  in  the  till  cell  of  the  grid  car¬ 
ries  with  it  the  weight 

u'm,(7i-  P i)~  “'„,,(  v, )  =  ( —  1 )”  exp(-y,)  f.„( 2 v, )  n/i  • 

Starling  with  ((/,,/>,)  and  the  initial  values  of  the  transla¬ 
tional  coordinates  Q  and  mmueiil.i  we  now  numcneally 


integrate  Hamilton’s  equations  of  motion  until  the  final 
condition  is  reached.  As  in  the  usual  QC  method,  we 
then  obtain  the  final  energy  of  the  oscillator  or  the  cor¬ 
responding  final  y  value  y,1/1. 

The  collection  of  these  final  v  values,  with  the  weight 
factor  (*)„„( y,)  now  associated  with  y,v\  may  be  considered 
to  represent  the  final  Wigner  distribution  H'(  y,  I  =«). 

The  transition  probability  n-  m  can  then  be  obtained  by 

Pm-n  =  2rrAf  f  dq  f  dp  wmm(q ,  p)  W(q,  p,t  —) 

“ K'f  dyumm(y)W{y,l‘">) 

w-(y.tf>)  ««.(>■,)  .  (2D 

i 

where  K  and  K'  are  appropriate  normalization  constants. 

We  note  here  that  the  final  Wigner  distribution  is  con¬ 
sidered  to  be  a  function  of  y,  not  of  q  and  p  separately. 

If  we  were  to  consider  q  and  p  as  two  independent  vari¬ 
ables,  the  physical  distribution  obtained  would  depend 
on  the  choice  of  the  final  time.  However,  since  the  final 
state  is  in  general  a  superposition  of  energy  eigenstates, 
the  correct  final  Wigner  distribution  function  should  con¬ 
tain  off-diagonal  elements  uu,  i*j,  and  therefore  should 
depend  on  q  and  p  separately.  In  other  words,  all  infor¬ 
mation  concerning  the  off-diagonal  elements  u.'(y,  i*  /  ,  in 
the  Wigner  distribution  is  lost  in  our  treatment.  This  is 
the  price  that  we  have  to  pay  for  using  classical  equations 
of  motion  to  describe  the  dynamics.  The  same  difficulty 
also  exists  in  the  initial  Wigner  distribution.  One  might 
argue  that  the  initial  distribution  represented  by  the  set 
of  points  on  the  grid  is  meaningless,  in  the  sense  that 
the  distribution  would  be  physically  different  for  a  dif¬ 
ferent  choice  of  the  initial  time.  This  difficulty,  how¬ 
ever,  can  be  avoided  at  least  numerically  by  taking  suf¬ 
ficiently  many  points  on  the  grid  such  that  the  difference 
is  negligible.  This  is  ensured  in  our  calculations  by  in¬ 
creasing  the  number  of  points  on  the  grid  until  the  final 
results  change  little  with  a  further  increase. 

Calculations  have  been  carried  out  for  transition  prob¬ 
abilities  of  the  He-Hj  system,  using  the  proposed  S QC 
method  as  well  as  the  usual  QC  method.  The  results 
are  shown  in  Table  I  and  Fig.  1  for  various  initial  trans¬ 
lational  and  vibrational  energies,  along  with  the  “exact” 
results  of  Secrest  and  Johnson. 13  The  SQC  method  is 
seen  to  yield  reasonably  accurate  probabilities  for  the 
case  n  =  0,  i.e.,  when  the  oscillator  is  initially  in  its 
ground  state.  However,  for  ti  =  2  the  SQC  method  is  al¬ 
most  as  “inaccurate”  as  the  QC  method.  It  is  not  im¬ 
mediately  clear  why  the  SQC  method  becomes  less  ac¬ 
curate  as  v  increases.  The  most  interesting  feature  in 
our  results  is  that  tiic  .SQC  method  yields  Unite  proba¬ 
bilities  for  the  two  classically  forbidden  processes; 

0-  0  for  the  case  K  -  20  (in  units  of  and  1-3  for 

E  -  12  (see  Table  I).  The  same  method,  however,  yields 
zero  probability,  within  the  accuracy  of  the  method,  for 
the  other  two  classically  forbidden  processes;  2—  5  for 
E  -  16  and  2-  4  for  E  -  12.  This  only  means  that  the 
SQC  method,  as  developed  here,  is  capable  of  describ¬ 
ing  classically  forbidden  processes,  but  needs  further 
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FIG.  1.  Transition  probabilities  P  as  functions  of  final  vibra¬ 
tional  quantum  numbers  m,  calculated  from  the  quantum  me¬ 
chanical  (solid  linesl  (Hef.  13),  quasiclassical  (dashed  tinesl, 
and  the  statistical  quasiclassical  (dotted  linesl  methods.  Ini¬ 
tial  translational  and  vibrational  energies  £,  and  £„  are  given 
in  units  of  jAw,.  where  u„  is  the  vibrational  frequency  of  lb. 


improvement  to  be  sufficiently  accurate.  In  order  to 
further  study  the  accuracy  of  the  SQC  method  in  de¬ 
scribing  classically  forbidden  processes,  we  also  made 
calculations  for  the  case  E  -  8.  In  this  case  the  initial 
energy  is  so  low  that  most  transitions  arc  classically  for¬ 
bidden.  The  results  of  the  SQC  and  QC  calculations  are 
shown  in  Table  II,  along  with  the  QM  results  of  Secrest 
and  Johnson. 13  One  sees  that  the  SQC  method  is  very 
accurate  in  this  case.  This  is  made  more  evident  in 
Table  III.  Here  the  results  of  SQC  and  QC  calculations 
are  shown  for  a  collinear  He-HBr  collision  with  £  =  12 
(in  units  of  »/iu>0,  where  u.'0  is  now  the  vibrational  fre¬ 
quency  of  HBr).  Also  shown  are  the  corresponding  QM 
data  of  Clark  and  Dickinson. 14  In  this  case  all  transi¬ 
tions  are  classically  forbidden,  as  can  be  seen  from  the 
QC  data.  Nevertheless,  the  SQC  method  gives  quite  ac¬ 
curate  probabilities  for  all  the  forbidden  processes  listed 
in  Table  III.  On  the  other  hand,  the  SQC  method  cannot 
accurately  describe  processes  that  have  probabilities 
less  than  ~  10‘3.  For  example,  the  QM  calculations" 
show  that  the  transition  0-2  (and  2-0)  has  the  proba¬ 
bility  ~  10".  This,  however,  is  beyond  the  accuracy  of 
the  SQC  method  in  its  present  form. 

A  useful  criterion  by  which  to  judge  the  accuracy  of  a 
method  is  whether  forward  and  reverse  probabilities 
computed  using  that  method  are  equal  at  the  same  total 
energy;  i.e.,  whether  the  transition  probabilities  obey 
microscopic  reversibility.  QM  probabilities  obey  this 
principle,  whereas  QC  ones  often  badly  violate  it.  as  can 
be  seen  in  Tables  I  and  II.  Tables  I,  II,  and  III  also 


TABI.F  I.  Transition  probabilities  for  a  collinear  lle-li;  collision  calculated  from  the 
quantum  mechanical  (QM)*  (Kef.  1.1),  quasiclassical  iqC),  and  the  proposed  statistical 
quasiclassical  (SQC)  methods.  The  total  initial  energy  £  is  measured  in  units  ol 
where  o>„  is  the  vibrational  frequency  of  II,. 


QM* 

QC* 

sqc* 

QMC 

t,*rc 

SQC* 

Q.\rd 

QC* 

SQC* 

0  —  0 

(0.060)* 

0* 

0.016* 

(0.204) 

0.  300 

0.  187 

(0.  538) 

0.500 

0.520 

1 

0.218 

0.375 

0.  202 

0.434 

0.  275 

0.422 

0.  394 

0.  425 

0.412 

2 

0.366 

0.200 

0.351 

0.291 

0.313 

0.314 

U.  008 

0.  07 5 

0.008 

3 

0.267 

0.250 

0.  294 

0.071 

0.113 

0.077 

4 

0.089 

0.  175 

0.  106 

1  —  0 

0.218 

0.  263 

0. 199 

0.431 

0.350 

0.420 

0.394 

0.388 

0.411 

1 

(0.  286) 

0. 103 

0.285 

(0.0341 

0.  175 

0.065 

(0.  211) 

0.225 

0.  176 

2 

0.009 

0.  125 

0.012 

0.  220 

0.  163 

0.  151 

0.345 

0. 388 

0.385 

3 

0. 170 

0.  125 

0.001 

0.  261 

0.  238 

0.302 

0.  037* 

0* 

0.  028* 

4 

0.210 

0.175 

0.202 

0.051 

0.075 

0.061 

5 

0.077 

0.  150 

0.  121 

2-0 

0.  366 

0.  275 

0.346 

0.291 

0.300 

0.309 

0.  008 

0.  100 

0.065 

1 

0.001) 

0.  150 

0.038 

0.220 

0.  175 

0.  163 

0.345 

0.  350 

0.371 

2 

(0.  207) 

0.  100 

0.  143 

(0.034) 

0.  125 

0.  101 

(0.34M 

0.  250 

0.  280 

3 

0.018 

0.  125 

0.  100 

0.  250 

0.  150 

0.  157 

0.  233 

0.300 

0.  281 

4 

0.  169 

0.  too 

0.052 

0.189 

0.  250 

0.270 

<1.  006* 

0* 

0* 

5 

0.  191 

0.200 

0.  245 

0.016* 

0* 

0* 

6 

0.037 

0.050 

0.  oos 

‘Quantum  mechanical  probabilities  fur  elastic  transitions  given  In  parentheses  are  up|>cr 
bounds. 

*A  j  20. 
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TABLE  II.  Trinnition  prolubdliies  for 
(  colllnear  He- II,  collision  calculntml 
from  the  quantum  mechanical  HfMi* 
(Ref.  13),  quasiclattHical  IQL),  and  the 
statistical  quasiclassical  ISQL)  method* 


1  1  ■'  ~s  = 

QM* 

QL* 

St*’* 

0  —  0 

(0.  *92) 

0.73H 

0.  893 

1 

0.  108 

0.263 

0.  107 

1-0 

0.  10* 

0.  238 

0.  106 

l 

(0.  *501 

0.888 

0.  863 

*> 

0.0*2 

0.075 

0.031 

2  -  1 

0.  0*2* 

0* 

0.033* 

■I 

(0.9571 

1 

0.  967 

3 

0. 001* 

0* 

0* 

*Quantui 

n\  mechanic 

al  probabilities  lor 

•llstic 

transitions 

given  in  \ 

•arenthese? 

are  up|*er  Itotmds . 

*£  ’  i  tn 

units  of  \fi 

a where 

,  is  the 

vibrational  Irequencs  ul  II  . 

Vl,3sh-,ll\  (ori-idilen  pn  kw-sm-*  . 

show  that  SQC  probabilities  do  obey  iiiicroMopic  re¬ 
versibility  to  .1  good  aii|)ro\ii:iation.  This  is  one  gootl 
indication  that  the  SQL  method  is  superior  tu  the  QC 
method. 

We  have  stiown  th.it  the  SQL  method  yields  more  ac¬ 
curate  probabilities  than  the  QL  method  lor  volliue.ir 
He-H,  and  lle-HHr  collisions .  One  -ti ould  Keep  in  mind, 
however,  that  the  QL  method  fares  nm.it  orablv  fur  Hie 
particular  systems  chosen  for  our  calculations.  Sys¬ 
tems  with  low  particle  masses  olfer  a  stringent  lest  of 
the  QC  method.  furthermore,  quantum  tunnel  me.  el  fects 
manifest  themselves  more  strongly  hi  a  collmear  colli - 
Sion  than  in  a  realistic  three-dimensional  collision,  1  or 
systems  in  w  hich  quantum  effects  are  small  and  the  QC 
method  fares  well,  ttie  use  of  the  SQL  method  may  not 
be  fully  justified.  This  is  pa rticula rly  s.i  because  the 
SQC  metliod  requires  a  l.irc.e  number  of  traioctories  — 
typically  an  order  ol  magnitude  more  than  required  hi 
the  QC  method — to  be  computed,  as  each  element  ) 

of  the  Winner  distribution  must  be  propagated.  In  fact, 
this  additional  computational  effort  represents  a  major 
disadvantage  of  the  SQC  method  as  compared  with  the 
QC  method.  However,  the  SQC  method  still  requires 
much  U:s s  com put.it ion.! I  eliort  ih.in  the  QM  method, 

©spec tally  when  the  collision  energy  is  high  and  a  large 
nu ml h? r  of  channels  need  to  l>e  considered.  This  is  he 
cause  the  number  of  trajectories  requ  red  to  be  com¬ 
puted  m  the  SQc  method  does  not  increase  rapidly  with 
respect  to  the  collision  energy,  whereas  the  QM  method 
goes  rapidly  out  ot  reach  as  the  collision  energy  is  in¬ 
creased.  Perhaps  the  strongest  point  m  favor  of  the 


TABLE  III.  TranalUon  probabilities  for  a 
col  linear  lle-IIBr  collision  calculated  from 
the  quantum  mechanical  KjM)*  (Hef,  1*0, 
quiiiclantilcai  and  the  statistical 

quaatclasalcal  <St}C)  methods. 


QM* 

QL* 

SQL* 

0-0 

(0.976) 

1 

0.975 

1 

0.  02** 

0* 

0.  025* 

1-0 

0.02** 

0* 

0.025* 

1 

(0.960) 

1 

0.956 

2 

0.016* 

0* 

0.019* 

2-1 

0.016* 

0* 

0.019* 

2 

10.9*0) 

1 

0.  97* 

3 

0.005* 

0* 

0.003* 

•Quantum  mechanical  probabilities  tor 
elastic  transitions  *;ivcn  in  parentheses 
are  up| wr  bounds. 

-  12  in  units  of  jfi u  0.  where  v„  is  the 
vibrational  frei|uenc.  of  MPr. 

°L'i assiL' a ! I \  forbidden  processes. 

SQC  metliod  is  that  quantum  tunneling  effects  are  de¬ 
scribed  with  real  classical  trajectories.  This  represents 
an  important  improvement  over  the  scmiclas.weal  colli¬ 
sion  theory5  m  which  one  carries  the  burden  of  working 
with  complex-valued  trajectories. 
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A  new  semiclassical  decoupling  scheme  for  electronic 
transitions  in  molecular  collisions:  Application  to  vibrational 
to-electronic  energy  transfer 

Hai-Woong  Lee 
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A  new  semiclassical  decoupling  scheme  (the  trajectory -based  decoupling  scheme)  is  introduced  in  a 
computational  study  of  \ibrat tonal •tivelcctromc  energy  transfer  for  a  simple  model  system  that  simulates 
collinear  atom -diatom  collisions  The  probability  of  energy  transfer  (/*)  is  calculated  ipiasiciassically 
using  the  new  scheme  as  well  as  quantum  mechanically  as  a  function  of  the  atomic  electronic -energy 
separation  (A),  with  overall  givxl  agreement  between  the  two  sets  of  results  Classical  mechanics  with  (he 
new  decoupling  scheme  is  found  to  be  capable  of  predicting  resonance  behavior  whereas  an  earlier 
decoupling  scheme  (the  coordinate-based  decoupling  scheme)  failed  Interference  effects  are  not  exhibited 
in  P  vs  A  results 


*.  INTRODUCTION 

Vibrational -electronic  energy  transfer  processes 
such  as 

A  <  BC(r)  :A’«  BC(  r)  (1) 

are  of  fundamental  importance  m  the  theory  of  molecu¬ 
lar  collision  dynamics  l>ocause  they  represent  simple 
examples  of  electronically  nonadiabatic  collision  pro¬ 
cesses. 1  An  accurate  description  of  Process  il),  how¬ 
ever,  has  not  been  without  difficulties.  The  quantum 
mechanical  (QM >  coupled -channel  approach,  although 
most  accurate  and  reliable  of  all  existing  methods,  suf¬ 
fers  from  excessive  computer  time  and  expenses. 
Seiiiiclassie.il  (SC)  and  quasirlassieal  (QCI  methods,  as 
they  have  been  further  developed  for  the  description  of 
nonadiabatic  processes, :  are  considered  to  give  reliable 
results  in  mam  cases,  although  not  as  accurate  as  the 
QM  method.  However,  as  thev  treat  internal  molecular 
(vibrational  and  rotational)  degrees  ol  freedom  classical¬ 
ly,  both  methods  are  generally  lielicved  to  be  incapable 
of  describing  near  •  resonance  processes  involving  vibra¬ 
tional  or  rotational  energy  levels.  For  example,  SC 
calculations’  attempted  on  the  near -resonance  collinear 
collision  process 

MiM’/’,  ,)  •  H:(>,  0)  -  |lu!/')/:]  •  ||,(i,  1)  (2) 

yielded  energy  transfer  probabilities  in  strong  disagree¬ 
ment  with  the  QM  results.’  Thus,  it  seems  that  one  has 
to  relv  only  on  fullv  QM  calculations  to  correctly  dr- 
srritic  near  -  resonance  processes.  In  view  of  this  rather 
unfortunate  situation,  further  development  oil  the  SC  and 
QC  methods  wrms  ilesir.ihlc 

,,||fM'Jll  II  . . I  >  it  I  *\  till  All  I  <  •!  i  •  <  ’ll  It  «  >  >t  St  1 1  (it  it  It  Hi 

•  i  if.  Ii  tindt  i  f.i  ml  N"  A  I  •  •  •  I  . 

*'AII»ivl  1’  •  iii-ti  I*  !!•■»% ,  I*»74i  s»ij  i  .  1 1 1 1  ■  I  It  ,m«i 

It*  nit.  I  M  •  \  !■  i  1 1  i.  In  i  N  Ii  -I  ir.  f»V  . 
f,||»  si-  ir  i  »i  -sifc-n  "it'll*',  thi  Air  I  or.  r  <  M I  |t«-  nl  N  lent  It  |.‘  t!«- 

M  .lit  Ii  nit* I*  a  i  "itt i  . it  t  N"  |  I'l'iJH  ?s  t  llit  Vitimi. il 

A«  t oii.i'itii  s  -nt*<  S|Mt  t  Attmini  .1 1  .il ion  umit  i  ( >  i  .nit  No  Nm  , 

ITJ*',  .tint  Hit  N  itioii.il  .  ,  I.  on  )  oumi.iliun  umli  i  lii.mlNo. 

C  HI. 77  2  •  HJi. . 


On  the  other  hand,  one  might  argue  that  it  is  not  en¬ 
tirely  clear  that  the  resonance  effects  should  be  absent 
only  because  internal  motions  are  treated  classically. 

The  disagreement  between  the  QM  and  SC  calculations  on 
Process  (2)  may  be  due  not  to  the  failure  inherent  in  the 
SC  method  but  to  the  “decoupling”  approximation3  made 
on  the  nonadiabatic  coupling.  In  other  words,  one  may- 
only  need  to  improve  upon  the  approximation  concerning 
the  nonadiabatic  coupling  and  still  be  able  to  use  the  SC 
or  QC  method  for  the  description  of  near-resonance  pro¬ 
cesses.  This  point  deserves  further  investigation  be¬ 
cause  the  SC  and  in  particular  the  QC  method  represent 
niaior  computational  simplifications  over  the  QM  method. 
Miller3  recently  argueu  that,  using  an  oversimplified 
model,  the  SC  method  (but  not  the  QC  method)  should  be 
able  to  correctly  describe  near- resonance  collision  pro¬ 
cesses.  Unfortunately,  no  calculations  have  been  car¬ 
ried  out  on  realistic  models  or  systems. 

In  this  paper,  we  present  a  computational  study  of  the 
resonance  effects  in  the  simplest  of  vibrational -to-elee- 
tromc  (  V-H)  energy  transfer  processes  :  collinear  atom- 
harmonic  oscillator  collision.  Specifically ,  we  are  in¬ 
terested  m  the  process 

A  *  l!C(n  l)-A*.HC(i>  0),  (3) 

where  the  diatom  ftC  is  now  approximated  as  a  harmonic 
oscillator.  The  energy  transfer  probability  /'  for  Pro¬ 
cess  (3)  is  calculated  both  quantum  mechanically  and 
quasiclassically  as  the  electronic  energy  separation  » 
of  the  atom  A  is  varied  around  a  fixed  value  of  the  vibra¬ 
tional  energy  separation  of  the  harmonie  oscillator  IK’. 
The  purpose  of  our  study  is  to  investigate  the  resonant 
process,  and  m  particular,  to  propose  (and  test)  a  new 
decoupling  scheme  for  the  calculation  of  CJC  probabilities. 
This  new  scheme  will  lie  termed  the  trajectory -based  de¬ 
coupling  srlieme,  lor  reasons  to  lie  discussed  m  her.  IV. 
Hy  investigating  the  QC  pruu.iliililics  thus  calculated,  it 
is  hoped  that  we  ran  approach  an  aiiswci  as  to  whether 
the  t^C  (and  SC)  method  may  lie  used  to  study  near-reso¬ 
nance  processes. 


TOO  J  Oii.ii  Ptiyv  /3( ) )  i  July  1980 
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In  Sec.  II  the  atom-harmonic  oscillator  model  system 
ia  described.  Section  III  follows  with  a  description  of 
the  QM  approach  and  results,  and  the  QC  counterpart  is 
described  in  Sec.  IV  along  with  the  newly  proposed  de¬ 
coupling  scheme.  Results  ot  this  calculation  and  a  com¬ 
parison  between  this  new  decoupling  scheme  and  the  old 
one*  is  presented  in  Sec.  V. 

II.  THE  MODEL 

We  consider  the  collincar  energy  transfer  Process  (3) 
between  an  atom  A  and  a  diatom  DC  approximated  as  a 
harmonic  oscillator.  For  simplicity,  the  atom  A  is  as¬ 
sumed  to  interact  only  with  the  nearer  of  D  and  C,  and 
the  interaction  strength  to  be  exponentially  decreasing 
with  respect  to  the  internuelear  distance.  The  potential 
surfaces  are  chosen  to  mimic  those  used  for  the  Br-H2 
system  in  Ref.  4.  The  two  adiabatic  potential  surfaces 
IV,  and  Wt  are  given  by 

Wx H«-[(W«-W„)2  +  4tf2|l/2!  2  .  (4a) 

Hii+UHn- +4H\z)U'-}  2  .  (4b) 

where 

=A,  exp(  -  a(A  -  r-p0)l  ,  3  mw*(r- r,,)2  ,  5  X  ,  (5a) 

H**  -  Aj  exp[  -  o(A  -  ±  r-  p„)|  +  i  »/u>J(r-  rj'  4  ]  X  ,  (5b) 

and 

/2 

<5c) 

Here,  K  is  the  translational  coordinate  measured  from 
the  c.m.  of  BC  to  A,  and  r  is  the  vibrational  coordinate. 
The  first  terms  in  //,,  and  UZ2  represent  the  exponentially 
decreasing  interaction  strength  between  A  and  the  nearer 
of  B  and  C,  and  the  second  terms  the  harmonic  oscilla¬ 
tor  potential  function.  The  potential  parameters,  chosen 
to  resemble  those  of  the  Br-H2  system,4  are  given  as 

A,  =  3x  10*5,  A,  =  5*  lO'5,  u  3,  p0  4  , 

<ut-0.  02,  r„  1.4,  mA-  1.5x10s,  wB  -  wc  -  2000, 

where  niA,  >irB,  mc  arc  respectively  masses  of  A,  B, 
and  C.  Thus,  the  reduced  mass  »i  of  the  harmonic  os- 


02 


t  0.t‘ 

o 


~  I 

R  (Bohr) 

FKi.  I.  The  (r.'liiillil  i-.ii'inv  :;iir(aei'H  W ,  .mil  W.  vs  the  tramda- 
tlon.il  coordinate  H  I'hc  vibrational  ciHirdmalc  Ik  fixed  at  1.4 
and  A  Ik  iihHuninl  In*  r<|li.il  In  O.Olti. 


FIG.  2.  The*  nonadiabatic 
coupling  vs  R,  "  R  -  i  r.  X  is 
assumed  to  be  equal  to  0.  OKI 


dilator  is  given  by  1000.  All  quantities  throughout  this 
paper  are  given  in  atomic  units  unless  otherwise  stated. 
For  example,  A,  and  A,  arc  given  in  hartrees  and  p0  and 
r0  in  bohrs.  X  is  the  electronic  energy  separation  of  the 
two  electronic  levels  of  the  atom  A: 

lim  (W,  -  W,)  =  X  . 

*-* 

One  may  assume  that  the  two  levels  are  separated  by  the 
spin-orbit  interaction  us  in  “Br.  ”  However,  as  noted 
in  the  Introduction,  X  is  considered  to  be  a  variable  in 
our  computational  study  and  therefore  is  varied  around 
the  vibrational  energy  separation,  R*:a  -  0.02. 

In  Fig.  1,  the  potential  energy  surfaces  IV,  and  H;  arc 
drawn  with  respect  to  the  translational  coordinate  H  for 
a  fixed  vibrational  coordinate,  r  -  »•„  =  1. 4,  for  the  case 
x  -  0.016.  It  can  be  immediately  seen  that  the  potential 
difference  IV2  -  IV,  is  a  monotonically  decreasing  func¬ 
tion  of  R,  and  therefore  no  avoided  crossing  occurs  in 
our  model  system. 


Having  defined  the  potential  energy  surfa-es,  we  now 
turn  to  the  nonadiabatic  coupling.  The  nonadiabatic  cou¬ 
pling  is  given  by  the  matrix  elements* 


1 5  ? r,m  r*"  A H - /(Wt  ~  w,')  - 

(6a) 

and 

dr  -Y,  £ri-r*  £*-.  <»*'*-"',)  , 

(6b) 

where 

^ii 1  Tjj  -  H„/[  ( IV,  -  //,i I2  4  H23l1/,  , 

(7a) 

r„  *-rtl«  »„/[(«', -//„? ,  n\t\ul  . 

(7l>) 

A  straightforward  calculation  yields 

'U  -  2<  -^X*(V-;x)  «/,*X-)<'*{|V 

-«/  •  t  x*)l/-T 

<I'X*|,/:{(<1I  4  (</  4  |-X*),/3|2  4  £.X*},/J  . 

(8) 

where 

V  (/Ij -A,)exp[  -  o(R >-p„)|  .  j  X  . 

(9) 

It  should  bo  noted  that  if*  and  it,  dc|»cnd  only  on  K,  ft 
-  j  r,  the  separation  between  A  and  the  nearer  of  1)  and 

C,  noton  Hind  r  separately.  In  Fig.  2,  Id*!  is  plotted 
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FIG.  3.  The  QM  and  QC  probabilities  vs  X  for  the  V-E  energy 
transfer  process  (3).  Collision  energy  =  «.  053.  The  solid  curve 
represents  the  QM  probabilities  and  the  circles  represent  the 
QC  probabilities.  The  dotted  curve  is  a  smooth  lit  to  the  cir¬ 
cles.  The  triangles  represent  the  QC  probabilities  obtained 
using  the  decoupling  scheme  of  Koinornicki  et  a! . 


FIG.  a.  The  QM  and  QC  probabilities  vs  1  f..i  the  V—E  one ruv 
transfer  process  (3).  Collision  energv  -0.03j.  |  he  solul 

curve  represents  the  QM  probabilities  and  the  cm  .es  n  pie.-  i 
the  QC  probabilities.  The  dotted  curve  is  a  sir.  lit  to  tin 
circles. 


against  for  the  case  1  =  0.  016.  The  nonadiabatic  cou¬ 
pling  is  seen  to  be  appreciable  over  the  interval  of  about 
2  bohr. 


results  was  obtained.  The  results  art  shown  in  Figs 
3-5,  in  which  the  energy  transfer  prubatiiliti.  .*>  ar- 
plotted  with  respect  to  '  lor  fixed  toliiaion  ci.oigies 
The  following  observations  and  interpretations  ran  tie 
made. 


III.  QUANTUM  MECHANICAL  APPROACH 

In  this  section,  we  report  the  results  of  QM  calcula¬ 
tions  for  the  V-E  energy  transfer  process  ' 3 1  performed 
on  the  model  system  described  in  the  previous  section. 
The  solution  to  this  problem  was  obtained  in  the  usual 
manner  by  solving  a  set  of  close-coupled  differential 
equations,  subject  to  scattering  boundary  conditions. 
Since  the  interaction  potentials  arc  modeled  as  exponen¬ 
tial  functions,  the  required  matrix  elements  over  the 
harmonic  oscillator  basis  can  be  performed  analytical¬ 
ly.7  The  I.ight-Walker  /{-matrix  method’  was  employed 
to  solve  the  close-coupled  equations,  and  the  step  size 
and  the  initial  and  final  limits  of  integration  of  the  trans¬ 
lational  coordinate  were  varied  until  convergence  of  the 


02 


KKi .  4.  The  QM  and  l,H  probabilities  vs  V  fur  llic  V-E  energv 
funster  process  (31.  Collision  rncrgv  n.ni.V  I  hr  solid 
curve  represent.*  llu-  QM  prob.ibi lilies  .uni  the  eireh-.s  repre 
Kent  the  QC  probabilities  'I  he  dotted  curve  is  a  smooth  lit  to 
the  circles. 


(1)  Each  curve  has  a  single  peak  structure  ref!,  rt. 
the  t'-£  resonance  effeet  under  consideration,  and  tt.< 
resonance  width  tends  to  •>  t  smaller  .it  lower  eoll'  -.oi 
energies  For  the  energy  r  in,  t  c  inside  red  ‘.ere,  tl,< 
width  is  of  the  order  of  0.  01  hartree  0.  27  «  \  . 

(2)  The  peaks  occur  at  values  of  *  smaller  than  the  \  i 
brutionul  energy  separatioii  0.02.  In  Table  l.  the  is  as 
positions  along  with  the  jx-as  probabilities  an  lisle  ,. 

Of  course,  one  would  not  ex|X>ct  the  peuits  to  occur  at 
exactlv  >  0.02  since  the  energy  levels  arc  perturbed 

by  the  collision.  In  lari,  the  rurves  are  strong,',  asv  r, 
metric  aliout  '  0.  u2;  the  probabilities  decrease  rap.nl 

as  x  is  increased  from  0.02  while  they  go  throu  li  a  i«-.i 
as  x  is  decreased  from  0  02 .  This  behavior  ran  I*  un¬ 
derstood  if  we  note  that  the  two  potential  surlares  rep¬ 
resenting  the  initial  and  llic  final  quantum  state:  do  not 
cross  If  >  is  greater  than  0,u2.  On  fhc  other  hand,  U  x 
is  smaller  than  0.02,  they  do  cross  and  th"  transition  i* 
more  likely  to  be  made.  In  jiartirular,  one  might  ex- 
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pact  that  the  transition  probability  is  the  largest  if  the 
crossing  occurs  at  a  point  of  the  maximum  nonadiabatir 
coupling.  A  straightforward  calculation  shows  that  this 
condition  is  satisfied  if  X  -  0.  012.  According  to  Table  1,  _ 

the  actual  peaks  occur  at  the  values  of  X  between  0.  012,  £ 

a  value  derived  from  the  above  expectation,  and  0.  020, 
a  value  expected  from  the  unperturbed  energy  level 
structure. 

We  have  also  calculated  probabilities  for  the  transi¬ 

tion  process 

A  ♦  BC(r,  =  1 )  —  A*  *  BC(  r,  11  (10) 

using  the  same  QM  method.  The  branching  ratio  of  the 
nonresonant  process  (10)  to  the  near -resonant  process 
(3)  for  J  particular  collision  energy  is  shown  in  Table 
11  for  various  values  of  A. 

IV  OUASICLASSICAL  APPROACH  DECOUPLING 
SCHEME 

In  discussing  nonadiabatir  transitions  in  quasictassicai 
terms,  one  usually  refers  to  the  Landau-Zener  type  of 
approximations  at  (avoided)  crossings  of  two  surfaces. 

The  two  surfaces  It',  and  Vf2  of  our  model,  however,  do 
not  possess  an  avoided  crossing  as  seen  in  Sec.  II.  They 
nevertheless  possess  complex  branch  points  at  which  the 
two  surfaces  cross  provided  that  they  are  allowed  to  take 
on  complex  values  and  at  which,  in  the  SC  theory,  non¬ 
adiabatir  transitions  are  allowed  to  take  place.  We 
therefore  compute  transition  probabilities  at  the  branch 
points  using  the  method  normally  employed  in  the  SC 
theory.  However,  the  phase  terms  rosi>onsibfo  for  quan¬ 
tum  interference  efforts  arc  totally  neglected  and,  for 
this  reason,  our  approach  is  still  termed  quasiclassical* 
following  the  usual  convention. 

The  branch  points  of  our  model  system  can  lie  easily 
found  by  solving  the  equation 

0.  (ID 

where  R and  rt  are  now  considered  to  assume  complex 

values.  Substituting  Fqs.  (4)  and  (5)  into  Eq.  (Ill,  we 

obtain 

{(Aj-A.loxpl  -a<K# r. -p0)l*  1 '!!.|  '2  (12) 

which  immediately  yields  the  following  equation  for  the 
branch  points, 

"  ,  (13a) 

u  P,  -  *  In-:  -----  ,  (13b) 

a  \-  A, 

r  —  cos'1  \  ,  (13c) 

o 

Even  with  tins  simple  .ui.ilytie  sliuiture  of  the  branch 
points,  a  very  large  ntimU  r  of  computational  efforts 
arc  still  needed  to  obtain  transition  probabilities,  if  the 
full  coinpli  \  valu'  d  tr.ip'clnrv  method  of  (he  SC  theory 
Is  to  lie  used.  A  considerable  computational  simplifica¬ 
tion  can  lie  achieved,  however,  with  the  recently  pro¬ 
posed  decoupling  scheme,  1  Since  we  will  use  an  im¬ 
proved  version  of  tins  scheme  in  our  calculations,  a 
brief  review  of  tin  scheme  is  in  order. 


FICJ.  6.  Collision  process  viewed  In  the  R-r  plane.  The 
branch  points  lorm  a  straight  line.  For  the  ease  A  =  u.01(i,  the 
equation  for  the  line  is  R  -  j  r  =  1.77.  The  solid  and  dotted 
curves  represent  a  typical  trajectory  on  the  way  in  and  out, 
respectively.  Ihe  choice  of  the  coordinates  ')  and  ;  according 
to  the  newly  proposed  decoupling  scheme  is  illustrated  at  the 
point  Q . 


The  decoupling  scheme  first  invokes  the  localization 
of  nonadiabatic  transitions  by  assuming  that  the  trajec¬ 
tories  remain  on  the  real  axis  except  when  they  hit  a 
point  that  coincides  with  the  real  part  of  the  branch 
points.  It  then  assumes  that  all  but  one  appropriately 
chosen  coordinate  p  and  corresponding  momentum  P, 
can  be  "frozen"  at  this  point,  and  calculates  the  transi¬ 
tion  probability  in  a  StUckelberg-like  manner  allowing 
only  the  “decoupled"  coordinate  n  and  momentum  Pn  to 
be  analytically  continued  into  the  complex  plane.  The 
essence  of  the  decoupling  scheme  is  therefore  to  reduce 
a  problem  with  multidegrees  of  freedom  to  one  with  a 
single  degree  of  freedom  analogous  to  an  atom -at am 
collision. 

We  now  wish  to  ronsider  our  atom-harmonic  oscilla¬ 
tor  collision  process  w  ithin  the  decoupling  picture. 

There  are  two  degrees  of  freedom,  translational  and  vi¬ 
brational,  in  tins  rollinear  process  as  represented  by  the 
coordinates  R  and  r  (and  corresponding  momenta  PR  and 
Pr).  For  this  case,  it  is  most  convenient  to  use  the 
R-r  plane  to  discuss  the  motion  of  the  system.  In  this 
plane,  the  real  parts  of  the  branch  points  given  by  Eq. 
(13b)  forms  a  straight  line,  as  shown  in  Fig.  6  for  the 
case  X  -  0.016.  The  system  moves  along  one  of  the  sur¬ 
faces,  or  It),  according  to  Hamilton's  equations  of 
motion,  until  it  encounters  a  point  on  the  line.  A  typical 
trajectory  is  shown  in  Fig.  G,  in  which  the  line  of  the 
branch  points  is  encountered  twice  during  the  eollision 
proress.  Whieh  points  on  the  line  the  trajectory  hap¬ 
pens  to  hit,  of  course,  depends  on  the  initial  conditions. 
To  make  our  discussion  definite,  let  us  consider  the 
point  (/  in  Fig.  f>  which  can  be  represented  by  (Ru,  rQ) 
in  the  R-r  plane.  To  obtain  Hie  transition  probability 
at  tins  |>oint  using  the  decoupling  scheme,  we  first  need 
to  rhoose  a  decoupled  coordinate  >)  and  a  frozen  coordi¬ 
nate  i,  each  of  which  can  lie  in  general  a  linear  combina¬ 
tion  of  It  and  ».  In  the  .scheme  promised  by  Komorneki 
rt  nl. ,  ’  i;  is  chosen  to  lie  It  and  (  fo  lie  r.  The  reason 
for  the  choiee  is  that  it  is  generally  believed  that  the 
component  of  the  momentum  most  effective  in  leading  to 
nonadiabatic  transitions  is  the  one  perpendicular  to  the 
line  of  tlu.'  real  parts  of  the  branch  |k>iii1s,  and  in  their 
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TABLE  II.  The  QM  and  QC  branching 
ratios  of  Process  (10)  to  Process  (3). 


X 

QM* 

QCb 

QCc 

0.012 

0.666 

1.23 

18.  7 

0.014 

0.405 

0.727 

5.23 

0.016 

0.266 

0.511 

4.82 

0.018 

0. 192 

0.463 

1.42 

0.020 

0. 156 

0.391 

5.63 

0.022 

0. 147 

0.254 

2.56 

•Collision  energy  =  0.  055. 

^Collision  energy  =  0. 055.  The  new  de¬ 
coupling  scheme  is  used. 

'Collision  energy  -  0.  055.  The  decou¬ 
pling  of  Komornicki  et  al .  is  used. 


study,  this  happens  to  be  very  nearly  the  translational 
momentum.  Within  this  scheme  the  transition  proba¬ 
bility  is  given  by 

p=ex p(-26)  ,  (14) 

where 


6=I*n  f*<,'‘'dR[fel(K)-t2(K)|  ,  (15) 

k,(R)  =  {2M[E-Wl(R,  rQ)l}1'2  ,  j  -1,  2  ,  (16) 

A/=  »)A(»uB  +  »ic)/(mA*  w,  +  x/c)  ,  (17) 

and  E  is  the  total  energy.  While  the  translational  coor¬ 
dinate  R  and  momentum  PK  are  allowed  to  proceed  in  the 
complex  plane,  the  vibrational  coordinate  r  and  momen¬ 
tum  P,  are  assumed  to  be  frozen  at  the  values  at  Q, 


r=rQ,  Pr  =  PTq  .  (18) 

Thus,  only  the  component  of  momentum  in  the  R  direc¬ 
tion  will  be  changed  after  the  transition  is  made, 


P\ 

2.1/ 


4  W, 


p'l 

W, +  H> 


(19) 


where  the  prime  is  used  to  denote  the  quantities  after 
the  transition.  The  above  choice  of  the  decoupled  and 
frozen  coordinates,  however,  presents  two  immediate 
difficulties.  First,  the  laws  of  classical  mechanics  are 
violated  by  assuming  that  both  coordinate  and  momen¬ 
tum  are  frozen  at  finite  nonaeio  values  as  dictated  by 
Eq.  (18).  The  position  r  can  be  fixed  at  rQ  only  if  Prq 
is  zero.  Second,  Eq.  (19)  suggests  that  the  change  in 
electronic  energy  represented  by  the  potential  jump 
W,-  Hj.  is  accompanied  by  the  change  in  translational 
energy  rather  than  vibrational  energy.  This  may  incor¬ 
rectly  suppress  V-E  energy  transfer  because  the  QM 
data  in  Table  II  show  that  V-F.  energy  transfer  lakes 
place  more  easily  than  Iranslational-to-electronic  ( T-E ) 
energy  transfer  in  near- resonance  processes.  It  now  is 
not  surprising  that  the  E—  V  energy  transfer  probabilities 
in  the  Ur— II.  system  as  calculated  by  Komornicki  et  al. 
were  anomalously  small.  It  ala)  is  evident  that  a  more 
reasonable  ■•heiee  of  the  coordinates  ?;  and  £  will  give  a 
bctlcr  result.  We  therefore  priqiose  the  following  in  an 
attempt  to  overcome  the  above  difficulties. 

We  choose  i)  to  lie  the  coordinate  along  the  direction 


of  motion  and  £  perpendicular  to  it.  This  choice  is  il¬ 
lustrated  in  Fig.  6  at  the  point  Q  for  the  trajectory 
shown.  Here,  n  and  £  take  the  form 

r?  =  fi-0Qr,  (20a) 

t*R+r/fi0  ,  (20b) 

where 

0o  =  -VKo  .  (21) 

and  r0  and  R0  are  the  vibrational  and  translational  ve¬ 
locities,  respectively,  at  Q.  We  note  that  rj  and  £  are 
in  general  different  for  different  trajectories  and  there¬ 
fore  have  to  be  defined  for  each  trajectory  at  each  en¬ 
counter  with  a  branch  point.  For  this  reason,  our  new 
decoupling  scheme  may  be  termed  the  trajectory-based 
decoupling  scheme,  as  opposed  to  that  of  Komornicki 
et  al. ,  which  is  based  on  the  choice  of  an  optimum  co¬ 
ordinate.  We  will  refer  to  the  earlier  scheme  as  the 
coordinate-based  decoupling  scheme. 

The  advantage  of  our  choice  of  ))  and  £  is  that  £  is  al¬ 
ways  zero  and  therefore  there  is  no  conceptual  difficulty- 
in  freezing  this  coordinate, 

4  =  40  =*o  + pt=°  •  (22) 

The  transition  probability  at  Q  is  given  by  Eq.  (14), 
where  6  now  becomes 

6=Imf°  dr1[kl (tj)  -  fr.,(7j)l  ,  (23) 

"o 

where 


*1 q  ~  Rq -0oro  >  (24) 

^(7)  =  {2p„0[£-tF/()?,  £„)]}‘'2  ,  >  =  1.2  ,  (25) 

the  reduced  mass  governing  the  motion  of  the  sys¬ 
tem  in  the  rj  direction,  is  shown  in  the  Appendix  to  be 
given  by 


„  AI/P*„*  m 

'  ( P0  4  1/^7  ’ 


(26) 


and  r,0  is  the  imaginary  part  of  the  branch  point  that 
must  satisfy  Eqs.  (13).  Simple  algebra  immediately 
yields 


1 


^  '  1 4  (ia/2 


(27) 


Equations  (20)— (27)  now  enable  us  to  calculate  for  tran¬ 
sition  probabilities  within  the  new  decoupling  scheme. 
Upon  the  completion  of  a  transition,  P,  is  given  a  new 
value  according  to  the  energy  conservation  relation 


(28) 


Because  we  are  interested  in  the  case  of  a  low  colli¬ 
sion  energy  (/:,  $  0.  055),  it  occasionally  happens  during 
our  QC  calculations  that  there  is  not  sufficient  energy 
available  for  an  upward  transition.  In  this  case,  the 
formula  introduced  by  l.aing  cl  al.9  is  employed  to  cal¬ 
culate  probabilities,  with  obvious  modifications  to  ac¬ 
count  for  the  fact  that  now  the  decoupled  coordinate  >| 
fs  the  one  to  be  integrated  over.  In  the  case  of  very 
high  collision  energies  such  that  vibrational  motion  is 
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relatively  weak,  we  expect  the  trajectory -based  de¬ 
coupling  scheme  to  approach  the  coordinate -based  one. 
Furthermore,  in  more  general  collision  conditions  not 
Involving  resonance  transfer,  the  two  decoupling  schemes 
may  also  be  equivalent. 

V.  RESULTS  AND  CONCLUSION 

The  method  used  for  the  QC  calculations  is  the  usual 
classical  trajectory  method10,11  modified  to  include  the 
calculation  of  nonadiabatic  transition  probabilities  with 
the  trajectory-based  decoupling  scheme.  The  Hamilton 
equations  of  motion  were  numerically  integrated  using 
the  fourth-order  Adams -Moulton  formulas  with  the 
fourth-order  Runge-Kutta  integration  being  employed  to 
initiate  the  procedure.  When  the  system  encounters  a 
branch  point,  the  transition  probability  was  computed 
as  described  in  the  previous  section.  The  complex  in¬ 
tegration  in  Eq.  (23)  was  carried  out  using  the  four- 
point  Gauss -Legendre  quadrature  formula.  Upon  com¬ 
pletion  of  the  probability  calculation,  the  integration  of 
Hamilton  equations  was  resumed.  For  a  path  starting  on 
a  new  surface,  the  Runge-Kutta  formula  was  again 
needed  to  initiate  the  integration  with  a  new  set  of  initial 
variables  before  the  Adams-Moulton  formula  could  be 
used.  The  system  typically  encountered  two  or  four 
branch  points  during  the  collision  process.  For  the 
case  of  four  branch  points,  eight  different  paths  must 
be  considered.  The  final  vibrational  quantum  number 
was  assigned  using  the  usual  quasiclassical  box.  The 
total  probability  for  the  energy  transfer  process  (3)  was 
then  obtained  from  the  equation 

J»=Z>r'Ar,  (29) 

i 

where  X  is  the  total  number  of  trajectories,  />,  is  the 
total  transition  probability  for  the  ith  trajectory  that  en¬ 
ters  the  box  v,-  0,  and  the  summation  runs  over  all  such 
trajectories. 

The  results  of  the  QC  calculations  are  shown  in  Figs. 

3,  4,  and  5.  The  circles  represent  the  QC  probabilities 
actually  obtained  from  the  calculations,  and  the  dotted 
curves  arc  smooth  fits  to  the  circles. 

In  Fig.  3,  the  result  of  QC  calculations  using  the  co- 
ordinate-based  decoupling  scheme  is  also  shown.  The 
probabilities,  represented  by  the  triangles,  are  rather 
unstable  with  respect  to  the  variation  of  A,  as  can  be 
seen  from  the  apparent  three-peak  structure.  Conse¬ 
quently,  it  is  hard  to  judge  whether  this  set  of  data  con¬ 
tains  any  resonance  effect  in  it.  The  superiority  of  the 
new  decoupling  scheme  can  be  more  clearly  seen  if  we 
also  consider  the  nonresonant  transition  (10).  In  Table 
II,  the  QC  branching  ratios  of  process  (10)  to  process 
(3)  obtained  with  the  two  decoupling  schemes  arc  shown 
along  with  the  QM  ratios.  The  QC  ratios  with  the  tra¬ 
jectory-based  scheme  show  a  reasonable  agreement  with 
the  QM  ratios.  However,  the  coordinate-based  scheme 
predicts  the  ratios  to  be  very  large  as  compared  to  the 
accurate  QM  ones.  This  seems  to  be  a  clear  indication 
tha't  the  coordinate-based  scheme  is  not  adequate  to  de¬ 
scribe  processes  involving  resonance  effects.  Figures 
3,  4,  and  5  also  exhibit  QM  results,  which  indicate  that 


the  QC  approach  with  the  trajectory-based  decoupling 
scheme  reproduces  the  resonance  peak  structure  rea¬ 
sonably  well. 

The  QC  peaks  occur  at  A  -  0. 012,  shifted  towards 
smaller  values  of  A  from  the  QM  peaks.  At  A  -  0.  012, 
the  potential  difference  tr,  -  if,  at  the  real  parts  of  the 
branch  points  is  equal  to  the  vibrational  energy  separa¬ 
tion  0.02;  i.e.,  this  A  fulfills  the  resonance  condition  at 
the  branch  points.  Speaking  in  QM  terms,  at  this  A,  the 
two  surfaces  representing  the  initial  and  final  quantum 
states  cross  at  the  real  parts  of  branch  points.  (See 
related  discussions  in  Sec.  III.  For  our  model,  those 
points  coincide  with  the  points  of  the  maximum  nonadi¬ 
abatic  coupling.)  Therefore,  it  is  reasonable  to  expect 
the  QC  peaks  to  be  at  A  -0.012  provided  that  the  reso¬ 
nance  effects  are  properly  taken  care  of.  Th...  the  QC 
peaks  occur  indeed  :<i  A  _  0.  012  indicates  that  the  reso¬ 
nance  effects  are  incorporated  in  our  QC  method  with  the 
new  decoupling  scheme.  The  accuracy  of  the  method, 
however,  is  limited  by  the  assumption  that  the  transition 
takes  place  only  at  branch  points.  This  assumption  is 
at  least  partly  responsible  for  the  inaccurate  position  ol 
the  QC  peaks.  Referring  to  Table  1,  we  see  that  the 
agreement  between  the  QC  and  the  QM  results  on  the 
position  of  the  peaks  is  better  at  higher  collision  ener¬ 
gies.  We  also  note  that  the  QC  probabilities  are  in  gen¬ 
eral  somewhat  larger  than  the  QM  probabilities.  One 
might  argue,  however,  that  the  inclusion  of  the  Niwitin 
factor9  will  bring  the  QC  probabilities  down  closer  to 
the  QM  ones. 

It  is  also  interesting  to  note  that  interference  effects, 
which  are  of  course  incorporated  in  the  QM  approarn 
but  suppressed  in  the  QC,  are  not  exhibited  in  the  transi¬ 
tion  probability  vs  A  results.  On  the  other  hand,  the 
absence  of  the  interference  effects  in  the  QC  method 
could  be  a  part  of  the  reason  why  the  QC  results  become 
less  accurate  at  lower  collision  energies.5 

More  stringent  tests  on  the  present  decoupling  scheme 
would  seem  to  rest  on  its  application  to  very  low  colli¬ 
sion  energies,  perhaps  even  to  cases  where  the  colli¬ 
sion  energy  is  much  lower  than  the  resonance  energy 
exchanged.  However,  it  is  debatable  whether  such  a 
condition  is  a  sine  qua  non  for  resonance  considerations. 
The  collision  energies  we  have  chosen  are,  in  fact, 
lower  than  or,  at  worst,  comparable  to  the  vibrational 
quantum  exchanged,  near  the  real  parts  of  the  branch 
point,  even  though  they  may  be  large  in  other  regions. 
Moreover,  our  exact  quantum  mechanical  results  also 
indicate  that  resonance  effects  do  persist  to  the  range 
of  collision  energies  considered  here,  even  though  these 
effects  tend  to  decrease  as  the  collision  energy  in¬ 
creases.  (The  single-peak  structure  of  the  resonant 
/’vs  A  curve  recurs  through  high  and  the  typical 
branching  ratio  for  resonant  to  nonresonant  transitions 
is  ~  10  at  the  peak  value  of  A,  for  fC  0.035.  )  Unfortu¬ 
nately  the  application  of  our  present  trajectory -based 
decoupling  scheme  is  restricted  to  cases  having  a  well- 
defined  real  velocity  at  the  real  parts  of  the  branch 
point,  which  is  not  likely  to  be  the  case  for  very  low 
collision  energies. 
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APPENDIX:  TRANSFORMATION  OF  COORDINATES 

In  terms  of  the  new  coordinates  r;  and  £  defined  by 
Eqs.  (20),  the  kinetic  energy  can  be  written  as12 

r=5  A/Rl  +  y  rnrz  =  \  p„)72  +  y  p,  k*  +  2  nvi  ■  (Al) 

Note  that  there  exists  a  cross  term  containing  t)£  be¬ 
cause  the  coordinates  rj  and  £  are  not  separable  in  gen¬ 
eral.  Substituting  Eqs.  (20)  into  Eq.  (A  1)  and  equating 
the  terms  that  have  the  same  powers  of  k  and  r,  we  ob¬ 
tain 


M  =  p,+  p, +  p  , 

(A2) 

m  =  p„|32  +  p,/|32-p  , 

(A3) 

0  =  -2m„P  +  2/x|/(3  +  h(1  >p-p)  . 

(A4) 

Solving  Eqs.  (A2)-(A4)  for  p,,  p,, 

and  ii,  we  obtain 

p„  =  (Af/02+m)/(i3+l'W  , 

(A  5) 

which  is  Eq.  (26)  in  Sec.  IV, 

4,  =  (;.$*+ w)/(i3+  1/UJ*  , 

(A6) 

and 

p  =  (2, If  -  2*w)/,(j3  +  1  '(3)2  . 

(A  7) 

Now  if  the  system  moves  under  the  constraint  that  £  re¬ 
mains  zero,  then 

T=|p„i72  =  P2/2p,,  ,  (AB) 

and 


*,=p,/A=(2p,r),/V/f=[2p,(£- w)Yn/n  ,  (ab> 

where  W  is  the  potential  energy. 


'For  a  review  of  the  theory  of  nonadiabatic  collision  processes, 
seej.  C.  Tulty,  in  Dynamics  of  Molecular  Collisions,  Modern 
Theoretical  Chemistry,  edited  by  VV.  H.  Miller  (Plenum,  New 
York,  1976),  Part  B,  p.  217. 

2Since  we  consider  nonadiabatic  processes  throughout  this 
paper,  we  will  be  using  the  term  “QC  method"  to  mean  the 
extension  of  the  usual  QC  method  to  account  tor  nonadiabatic 
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Pulsed  Visible  Chemiluminescence  Experiments 
M.  D.  Burrows  and  D.  J.  Benard 

Introduction 

In  order  to  achieve  lasing  on  an  electronic  transition  in  a  typical  dia¬ 
tomic  molecule  it  is  necessary  to  rapidly  generate  a  copious  supply  of 
molecules  in  the  upper  level  of  the  laser  transition.  At  shorter  wave¬ 
lengths  this  requirement  becomes  more  difficult  to  satisfy  because  of 
rapidly  decreasing  stimulated  emission  cross-sections.  Further,  when  a 
chemical  reaction  scheme  is  employed  as  the  source  of  excitation  it  is 
also  necessary  to  hold  down  the  reagent  concentrations  in  order  to  prevent 
quenching  of  the  lasing  species  by  the  reagent  or  product  molecules. 
Therefore,  it  is  requisite  in  such  cases  to  carry  out  the  chemical  reac¬ 
tions  in  as  rapid  a  way  as  possible  in  order  to  obtain  the  necessary  pump 
rate. 

Traditionally,  chemical  lasers  have  been  of  the  mixing  type  and  have 
employed  reactions  that  initiate  spontaneously  as  mixing  occurs.  This 
approach  is  most  appropriate  to  the  generation  of  CW  lasers,  however, 
since  the  mixing  is  limited  by  diffusion  which  often  proceeds  slowly 
compared  to  gas  kinetic  reaction  rates,  it  may  be  difficult  to  achieve 
lasing  on  diatomics  at  shorter  wavelengths  using  this  scheme.  As  an 
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alternative,  we  have  considered  laser  schemes  based  on  gas  phase  chemi¬ 
luminescence  reactions  that  do  not  initiate  spontaneously  but  proceed  via 
a  branched  chain  reaction  kinetics.  This  feature  not  only  allows  pre- 
mixing  of  the  reagents  but  also  provides  a  chemical  amplifier  effect  which 
allows  one  to  "trigger”  the  reaction  by  a  relatively  small  input  of  opti¬ 
cal  or  electrical  energy  which  produces  the  initial  chain  carrier  species. 
The  resultant  laser  action  would  therefore  be  pulsed,  which  may  be  advan¬ 
tageous  in  itself.  Enhanced  propagation  and  target  effects  have  been 
demonstrated  for  repetitively  pulsed  DF  chemical  lasers  relative  to  their 
CW  counterparts,'  however,  it  remains  to  be  seen  what  significant  advan¬ 
tages  may  accrue  for  pulsed  visible  lasers. 

Experimental  Method 

Fig.  1  shows  the  experimental  apparatus.  Chlorine  dioxide  (C102)  is 

generated  by  passing  Cl2  gas  over  NaC102  crystals.  (Ref.  1)  The  C102  is 

collected  in  a  cold  trap  prior  to  loading  into  the  reaction  cell.  A  fuel 

species,  typically  acetylene  (C2H2)  is  then  added  to  the  reaction  cell  and 

is  allowed  to  mix  with  the  chlorine  dioxide.  A  single  10  mJ,  20  nsec 

pulse  of  3500  A  light  from  a  XeF  excimer  laser  is  then  used  to  detonate 

the  mixture.  The  absorption  coefficient  of  C102  at  3500  X  is  so  high 

(3000  1/niole-cm)  that  the  10  mJ  pulse  is  almost  totally  absorbed  after 

traversing  less  than  1  cm  of  C102  at  partial  pressures  on  the  order  of  50 

Torr.  Irradiation  of  C102  at  this  wavelength  yields  predissociation  into 

C1CH0,  which  is  probably  followed  by  the  rapid  (10  1 1  cma/mol ecule-r.ec) 

* 

exothermic  reaction  CIO  +  0  -*  Cl  +  02.  (Ref.  2)  The  initial  UV  laser 
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Fig.  1.  Experimental  apparatus  for  pulsed  chemiluminescence  experi¬ 
ments.  The  excimcr  laser  can  be  operated  at  3500  X  (XcF) , 
2940  X  (KrF) ,  or  1950  X  (ArF) .  The  detector  consists  of 
either  an  interference  filter  and  photomultiplier  tube  or 
a  grating  monochromator  and  optical  multichannel  analyzer 
system.  The  quantum  efficiency  of  the  optical  multichannel 
analyzer  is  101,  decreasing  to  zero  at  the  sensitivity  limits 
of  3500  X  and  8400  X.  'Hie  resolution  is  5  X  per  channel  with 
a  memory  capacity  of  860  counts  for  each  channel. 


pulse,  therefore,  produces  high  densities  of  CIO,  0,  Cl,  and  excited  02  in 
a  small  volume  which  then  chemically  ignites  the  remaining  C102  and  fuel 
(C2H2).  The  detector  in  Fig.  1  consists  of  either  a  photomultiplier  tube 
and  interference  filter  for  intensity  versus  time  profiles  or  a  grating 
monochromator  and  optical  multichannel  analyzer  system  to  record  single 
detonation  spectra.  The  observation  side-arm  is  deliberately  located  a  . 
distance  away  from  the  end  window  so  that  optical  pumping  by  the  XeF  laser 
is  not  observed.  Typical  partial  pressures  are  150  Torr  of  C102  and  75  to 
100  Torr  of  C2H2.  Soot  formation  on  the  cell  windows  becomes  a  problem  if 
higher  C2H2  pressures  are  used. 

Results 

The  detonation  of  C102  and  C2II2  mixtures  resulted  in  an  intense  blue-green 
chemiluminescence.  A  typical  spectrum  of  the  chemiluminescence  is  shown 
in  Fig.  2.  The  observed  chemluminescence  is  easily  identified  as  the 
C2(d37i  *»  a3n)  Swan  bands.  (Ref.  3,  A)  The  emission  is  very  intense  and 
requires  the  use  of  neutral  density  filters  in  order  to  prevent  saturation 
of  the  optical  multichannel  analyzer  detector. 

Excitation  of  the  Swan  system  can  also  be  obtained  by  detonation  of 
HN3+C2H2+N20+He  mixtures  using  a  single  30  mJ,  20  nsec  pulse  of  2A90  A  KrF 
excimer  laser  light.  As  Fig.  3  shows,  strong  CN  Violet  band  emission  is 
also  observed  from  this  mixture.  Excitation  of  the  Swan  band  system  by 
such  mixtures  depends  much  more  critically  on  relative  reactant  concentra¬ 
tions,  but  produces  a  much  "cleaner"  and  better  developed  C2  Swan  system 
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Chemiluminescence  spectrum  from  the  detonation  of  a  mixture 
of  150  Torr  C107  plus  90  Torr  C7H, .  The  violet  band  heads 
at  5635,  5165,  and  4737  A  are  ttic‘"C2  Swan  bands  (d^ir->a-'n) . 
The  detonation  was  initiated  using  a  single  3500  A  excimer 
laser  pulse  and  the  spectrum  acquired  using  the  optical 
multichannel  analyzer  system  and  a  neutral  density  filter 
of  3. 
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Chemiluminescence  spectrum  from  the  detonation  of  a  mixture 
of  150  Torr  I1N3  +  50  Torr  C2I12  +  50  Torr  N20  +  200  Torr  He 
acquired  using  the  optical  multichannel  analyzer  system  and  a 
neutral  density  filter  of  2.  The  intense  violet  degraded  hand 
heads  at  5035,  5165,  and  4737  ft  are  the  C2  Swan  hands,  while 
the  violet  degraded  band  heads  at  4578,  4216,  and  3883  ft 
the  CN  violet  hands  .  The  weak  violet  degraded  head 

at  4314  ft  is  due  to  the  Cl  I  radical  A~A->X-ir  transition. 


than  does  the  detonation  of  C102  and  C2H2  mixtures,  where  the  Swan  bands 
•re  superimposed  on  a  strong  background  of  continuum  emission.  Detonation 
of  mixtures  containing  only  HN3  and  C2H2  produces  weak  Swan  bands,  while 
the  addition  of  He  buffer  gas  enhances  the  Swan  band  emission.  The  Swan 
band  intensity  is  increased  further  by  using  both  He  and  N20.  The  deton¬ 
ation  chemistry  of  HN3  and  C2H2  mixtures  appears  to  be  strongly  dependent 
on  trace  impurities  in  the  He  buffer  gas,  since  replacing  He  with  Ar 
completely  removes  both  the  CN  and  C2  emission.  Previous  researchers  have 
found  using  isotopic  substitution  that  the  carbon-carbon  bond  is  defin¬ 
itely  broken  in  reactions  producing  C2,  (Refs.  5,  6)  which  is  consistent 
with  the  observation  in  the  present  experiments  of  excited  CH  radicals 
(A2A  -*  X2/i  transition  at  4314  A)  when  mixtures  of  HN3+C2H2+N20+He  are 
detonated  by  KrF  laser  light. 

The  temporal  behavior  of  the  chemiluminescence  pulses  from  detonation  of 
C102  and  C2H2  mixtures  was  monitored  by  replacing  the  OMA  with  a  photo¬ 
multiplier  tube  and  5100  A  interference  filter.  A  typical  time  history  is 
shown  in  Fig.  4.  The  risetime  of  the  pulse  is  typically  30-50  usee,  while 
the  "burn-time"  exceeds  1  msec.  The  time  delay  before  the  rapid  onset  of 
light  emission  is  on  the  order  of  the  time  required  for  a  shock  front  to 
travel  from  the  laser  input  window  to  the  side-arm  where  light  emission  is 
observed . 


29 


5  mV /cm 


500  /ts*c/cm  200  fisec/cm 


ig.  4.  Time  history  of  emission  from  detonation  of  mixtures  of  50 

Torr  CIO?  and  50  Torr  C2II2  using  a  5100  A  interference  filter 
and  photomultiplier  tube. 


Discussion 


The  C2  Swan  bands  are  commonly  seen  in  hydrocarbon  flames  at  atmospheric 
pressure.  (Ref.  7).  The  d3rt  •+  a3n  transitions  are  highly  vertical  and 
thus  appear  in  a  characterstic  Av=0  vibrational  sequence  with  the  (0,0) 
band  head  at  5165  X.  Contributions  from  higher  vibrational  levels  appear 
at  shorter  wavelengths  spaced  by  u)^  ~  1800  cm  1 .  (Ref.  4)  Experiments  with 
C2H2  in  flames  have  shown  that  C2  Swan  band  emission  is  many  orders  of 
magnitude  more  intense  than  can  be  predicted  on  the  basis  of  the  adiabatic 
flame  temperature  and  statistical  (Boltzman)  population  distributions 
(Ref.  8)  Therefore  one  may  infer  that  the  C2  (d3n)  state  is  specifically 
populated  by  chemical  reactions  in  flames.  Further,  cavity  resonance 
experiments  have  failed  to  show  any  measureable  absorption  in  moderate 
pressure  flames  (Ref.  9)  from  which  one  might  surmise  that  C2(a37t)  is 
removed  via  a  chemical  process  since  statistically  the  C2(a37l)  state 
should  be  thermally  populated.  The  radiative  lifetime  of  the  C2(a3n) 
state  is  on  the  order  of  200  nsec.  (Ref.  10)  so  it  seems  probable  that 
rapid  chemical  removal  of  C2(a3n)  at  our  1019/cm3  concentrations  could 
sustain  an  inversion.  These  results  taken  in  conjunction  with  our 
experimental  observations  support  the  possibility  of  producing  a  pulsed 
visible  chemical  laser  on  the  C2  Swan  bands  in  premixed  C102  and  C2II2 
mixtures . 

At  this  Lime  the  mechanism  of  reaction  is  unknown,  thus  it  is  impossible 
to  speculate  in  any  detail  regarding  the  kinetics  of  the  reaction.  The 
simplest  interpretation  of  Fig.  4  (time  history  of  chemiluminescence)  is 


that  the  initial  rise  and  the  associated  high  frequency  structure  results 
from  a  chemically  driven  shock  wave  that  is  followed  by  uniform  burning  of 
the  reagents  behind  the  shock  front  which  persist  until  the  reactants  are 
depleted.  Thus  the  chemiluminescence  is  actually  a  quasi  CW  process  since 
the  reaction  time  is  long  compared  to  the  radiative  decay  rate.  One  would 
estimate  that  a  limiting  rate  of  ~  10  13  cm3/molecule-sec  controls  the 
chemiluminescence  process  based  on  the  reagent  concentration  and  the 
chemiluminescence  decay  rate. 

Of  course  whether  a  chemically  pumped  C2  laser  is  possible  will  depend 
critically  on  the  details  of  the  kinetics,  however,  one  may  show  that  if 
the  kinetics  efficiently  produce  C2(d37t)  and  selectively  destroy  C2(a3n) 
that  gain  coefficients  on  the  order  of  1/cm  can  be  easily  generated  by 
reacting  a  200  Torr  mixture  of  C102  and  C2H2  in  10  3  sec.  At  such  den¬ 
sities  the  chemical  energy  is  stored  at  the  rate  of  ~  1  KJ/liter  assuming 
efficient  production  and  extraction  of  the  chemical  energy  as  laser 
photons.  A  theoretical  specific  energy  in  excess  of  lMJ/lb  is  calculated 
by  dividing  the  photon  energy  by  the  mass  of  the  reagent  species.  In 
practice  it  is  unlikely,  however,  that  such  optimistic  figures  will 
actually  be  obtained;  however  as  a  result  of  the  thermal  expansion 
attendant  with  the  reactive  heat  release  it  would  appear  that  atmospheric 
pressure  recovery  is  feasible  in  a  ramjet  configuration. 

Summary 


We  have  shown  that  by  premixing  reagents  and  pulsing  a  chemiluminescent 
reaction  that  one  may  overcome  some  of  the  difficulties  inherent  to  the 
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production  of  a  visible  chemical  laser  operable  on  a  diatomic  species. 

The  reaction  of  CIO2  and  C2H2  is  one  such  reaction  that  produces  electron¬ 
ically  excited  C2  molecules  with  properties  appropriate  to  laser  action  on 
the  dsn  •*  asn  transitions,  given  kinetics  that  favor  inverson.  Such  lasers 
in  a  pulsed  mode  arc  capable  of  achieving  high  gain  coefficients  and  are 
characterized  by  high  specific  energy  and  high  energy  storage  per  unit 
volume  in  the  laser  cavity. 
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